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The Nambu-Bethe-Salpeter (NBS) wave function for two nucleons on the lattice has been shown 
to yield a non-local and energy-independent nucleon-nucleon (NN) potential, i/(r,r'). In prac- 
tice, the derivative expansion of t/(r,r') is currently employed to determine the potential at low 
energies. In this report, we study the magnitude of non-locality to check the convergence of such 
a derivative expansion. With quenched lattice QCD at m K = 530MeV, we compare the NN po- 
tentials at the center of mass energy E ~ MeV and at E ~ 45 MeV. We also investigate the 
angular momentum dependence of the spin singlet potential, by comparing the potentials in l So 
and D2 channels. We find that the non-locality and the angular momentum dependence in the 
above energy range are negligible within statistical errors. 
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1. Introduction 

The nucleon-nucleon (NN) potential, which describes the interaction between two nucleons 
at low energies, is widely used as a basic tool to study the structures and reactions of atomic 
nuclei. In the past few decades, several realistic NN potentials have been proposed: These are 
phenomenologically determined so as to reproduce the NN scattering phases of more than 4000 
data points with % 2 /N<\ i ~ 1 [0,^. §]■ Recently a first-principle method to extract the NN potential 
from the Nambu-Befhe-Salpeter (NBS) wave function in lattice QCD has been proposed by three 
of the present authors [Q,|5|]. This method can be generalized also to derive the hyperon potentials 
(YN and YY) [||, ^] and the three-nucleon potentials [||, which do not have enough experimental 
data. The potentials based on lattice QCD, once obtained in good accuracy, would become first- 
principle inputs to the studies of ordinary nuclei, hyper nuclei, and the interior of the neutron stars 
by using advanced techniques in nuclear many-body problems. 

The potential obtained from the NBS wave function in lattice QCD simulations is non-local 
and energy-independent, U (r,r'). To determine the potential at low energies, we have proposed to 
make a derivative expansion of U(r,r') and calculate the local coefficient functions successively 
[§]. The aim of this report is to examine the convergence of such a derivative expansion. In 
sec.§, we give a brief review of our method to construct the non-local potential and its derivative 
expansion. In sec.||, we show that the non-locality of U(f,f / ) can be checked by the center of 
mass energy (E) dependence and the orbital angular momentum (L) dependence of a "truncated 
local potential", V tmn (r;E,L). In sec.|], we show our numerical results of V tma (r;E,L = 0) at two 
different scattering energies, E ~ MeV and 45 MeV. In addition we examine the L dependence 
of the truncated local potential in the two spin-singlet channels, 'So and l Dj. We find negligible 
energy and angular momentum dependence of V tnm (r;E,L), which implies that the non-locality 
U (r, ?) is small at low energies up to E ~ 45 MeV. 

2. NN potential from Lattice QCD and its derivative expansion 

The equal-time NBS wave function 1 is defined as 

W> k ) = ^X>K* + rM*)l# = 2;*}, (2.1) 

X 

where "k" is the "asymptotic momentum" related to the total energy as W = 2^J rrv^ + k 2 . Also, 
the local composite operators for the nucleon are defined as np(y) = £ a bc{u a {y)Cy$db{y))d c p{y) 
and p a (x) = e a bc{u a (x)Cysdb{x))u ca {x), where a,b,c denote color indices, and C is the charge 
conjugation matrix. It has been discussed that <p(r,k) below the pion production threshold satisfies 
the Schrodinger-type equation [|],[5[]: 

(A + k 2 ) <p(r;k)=2n J d 3 r U(r,r) <j>(f ,k), (2.2) 

which defines the non-local potential and /^-independent potential U(f,r') with /I = being 
the reduced mass. By using the Nishijima-Zimmermann-Haag reduction formula for composite 

'The NBS wave function was also utilized to investigate Luscher's finite volume method for the pion-pion scattering 
phase shift in Ref . |TTJ] . 
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operators, it can be shown that, for large \r\, (j)(r;k) has an asymptotic form characterized by the 
scattering phase 8 (k) [||] : 

sm(kr + 8(k)) 

<Kr;k) { K JJ +■••■ (2.3) 

kr 

This functional form is exactly the same as the asymptotic form of the scattering wave in the 
non-relativistic quantum mechanics. Therefore our non-local potential U (r,r ! ) is to reproduce the 
scattering phase 8 (k) below the pion production threshold. 

Although the non-local potential 11(7,7') has nice formal properties, its explicit construction 
is not easy in lattice QCD simulations, where the number of available NBS wave functions for 
different scattering energies are limited. Therefore, we employ the derivative expansion as shown 
below and determine the local coefficient functions successively by using low-energy NBS wave 
functions: 



U(r,f , [ 



VZ(r) + V£(r) (a l -a 2 )+Vl(r) S n +V[ s (r)L-S+^(V 2 ) 



8(r-r'). (2.4) 



LO NLO 



Here / = 0, 1 denotes the total isospin of the two nucleons and S = (81 + 02) /2 is the total spin 
operator. Also, the tensor operator is written as S12 = 3(di • r) (02 • r)/r 2 — 0\ ■ The LO terms 
such as the central potential in the 1 Sq channel, and the central and the tensor potentials in the 
3 5i — 3 D\ channel have been obtained so far at CM energy E ~ MeV (see e.g. a recent summary 



3. Energy and angular momentum dependence as a test of non-locality 

In this section we discuss that the convergence of the derivative expansion of U(r,r') in the 
NN system can be examined by studying the energy and the angular momentum dependence of the 
truncated local potential at finite E. 2 For simplicity, we restrict ourselves to the spin-singlet case 
(S = 0) where Sn and L ■ S do not contribute. Then we obtain 

(A + k 2 )$ s =°(r;k)= 2,u [V c (r) + {V 2 , V p2 (r)} + V L ^r) L 2 + • • •] $ s=0 (r;k), (3.1) 

where Vc(r) = Vo(r) — 3V (r) is the LO potential in the derivative expansion, while the higher order 
terms with L 2n , V 2 " are referred to N 2,, LO terms. If we have independent NBS wave functions for 
different energies and angular momenta, we can determine N 2 "LO terms order by order starting 
from a small n. For such a procedure making sense, the magnitude of the potential at the leading 
order should be dominant at low energies and low angular momenta. A simplest way to check this 
is to define a "truncated local potential" given below and study its E and L dependence: 

If V^ an (r;E ,L) has small E and L dependence, it is a good approximation of Vc(r). On the other 
hand, if it has large E and/or L dependence, one cannot neglect V p 2 (r), V L i (r) etc, so that they must 
be determined together with Vfc(r) by using the NBS wave functions for several different values of 
E and L. 



2 In the Ising field theory, it is analytically shown that the energy-dependence is weak at low energy, indicating that 
the non-locality of the potential is weak WM. 
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4. Numerical Simulations and results 



4.1 Lattice QCD setup 

We employ the standard plaquette gauge action on a 32 3 x 48 lattice at j8 = 5.7 for quenched 
gauge configurations. Quark propagators are calculated with the standard Wilson quark action at 
fc = 0.1665. This setup gives the lattice spacing a -1 = 1.44(2) GeV (a ~ 0.137 fm) from m p input, 



the spatial extension L = 32a ~ 4.4 fm, m n ~ 0.53 GeV and ~ 1.33 GeV [13]. Quenched 
gauge configurations are generated by the heat bath algorithm with over relaxations. Potentials are 
measured with configurations separated by 200 sweep, and 4000 configurations are accumulated to 
obtain results in this report. These calculations are performed on Blue Gene/L at KEK. 

The NBS wave function is obtained from the nucleon four-point correlator in the large t region, 

G^(x,y,t,t ) = (0\n(y,t)p(x,t)J pn (t )\0) = Y,M0\n(y)p(x)\B = 2;W n )e- w ^-^ 

n 

-> A (0\n(y)p(x)\B = 2;W )e-W^>\ (4.1) 



where W„ = 2Jmjf + k% and A n = (B = 2;W„\J f pn \0). Here J pn {t§) = P(to)N(t Q ) is a source 
operator located at t = t with P a = e a ,b,c {UaCy 5 D b ) U ca and Np = E a ,b,c {U a Cy 5 D h ) D c p. Up and 
down quark operators associated with a source function f(x) specified later are denoted as U (t) = 
Lx M (/ >■*)/(-?) and D(t) = ^xd(t,x)f(x). By examining the t dependence of potentials, we confirm 
that the ground state saturation for potentials is achieved at t — to > 9. 



4.2 Details of the calculation 

The NN potentials at CM energy E = k 2 /(2ju) ~ MeV with \i = m N /2 is constructed by 
using the periodic boundary condition (PBC), which is imposed on the quark fields along spatial 
directions. Since spatial momentum of each quark is discretized as p\ ~ 2%nijh («, € Z) in the PBC, 
the ground state corresponds to pj ~ 0. For the PBC, we employ the wall source, i.e., f(r) = 1 with 
the Coulomb gauge fixing only at t = to to enhance the overlap with the ground state. The wall 
source provides us with the orbital contribution, by which we can only study S-wave (L = 0) 
for the spin singlet case with possible contaminations from higher partial waves with L > 4. To 
improve the statistics, we calculate the NBS wave function four times on each configuration by 
putting the source point at four different time slices. 

The NN potentials at CM energy E ^ is constructed by using the anti-periodic boundary 
condition (APBC), which is imposed on the quark fields along spatial directions. Since the spatial 
momentum of each quark is discretized as ~ (2n + \ )n/L («,• € Z) in the APBC, the ground 
state corresponds to pi ~ ±7t/L, whose energy is given by E ~ 3 x (7r/L) 2 /(2jU) ^ 45 MeV for 
L ~ 4.4 fm. For the APBC, we employ four types of "momentum-wall" sources, i.e., f(r) = 
cos((x + y + z)7l/L),cos((—x + y + z)7t/L),cos((x — y + z)7l/L), and cos((— x — y + z)n/L) with 
the Coulomb gauge fixing only at to to enhance the overlap with the ground state with pi ~ ±tc/L. 
In addition to the orbital contribution, our momentum-wall source provides us with the T 2 + 
orbital contribution, which makes it possible to study D-wave (L = 2) for spin singlet case with 
possible contaminations from higher partial waves with L > 4. 
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The NBS wave functions for the S-wave and the D-wave are extracted by projection operators, 

23 

Hnk-r) = £ I z (r) (*«■)* <t> s=0 (R-nk), (4-2) 



24 ■ 



where %^ denotes the character of the representation Y for the cubic group O (See Table.[l]), 7?, 
denotes one of the 24 elements of the cubic group, and dy denotes the dimension of the represen- 
tation, i.e., d A \ = 1 and dj2 = 3. The calculation performed with Y = A\ in the PBC provides us 
with the NBS wave function for 1 Sq at E ~ MeV, while the calculation with Y = A\ and r 2 in the 
APBC provides us with NBS wave functions for ^0 and ! D 2 at E ~ 45 MeV, respectively. These 
NBS wave functions are inserted into eq.(3.2) to obtain V^" n (r,E ,L), where k 2 /2/x in this report is 
estimated from the non-interacting nucleons on the lattice, i.e., k 2 /2/j. = MeV for the PBC, and 
k 2 /2pL = 3(n/L) 2 /2n = 45 MeV for the APBC. 



4.3 Results 

In the top panel of Fig. [T[ we compare V<? un (r;£,L = 0) in the l S Q channel obtained at E ~ 45 
MeV (red) with that at E ~ MeV (blue). All data are taken at t — to = 9, where the ground 
state saturation is numerically confirmed as mentioned before. We observe that the agreement of 
the potentials between the two energies is quite good within statistical errors. The comparison of 
V' ( f un (r;£',L) in the l S (L = 0) channel (red) with that in the 1 D 2 (L = 2) channel (cyan) obtained 
at E ~ 45 MeV is shown in the bottom panel of Fig. [|. We observe that the angular momentum 
dependence of V^ un (r;E,L) is weak, though the statistical errors are rather large in the l Dj channel. 

By these comparisons, we find that the energy and the angular momentum dependence of 
VQ Un (r;E,L) is very weak within statistical errors. We therefore conclude that contributions from 
higher order terms in the derivative expansion are small and that the LO local potential in the 
derivative expansion is the good approximation for the non-local potential U (r,r') at least up to the 
energy E ~ 45 MeV and angular momentum L < 2. 



5. Summary and discussion 

We have studied the validity of the derivative expansion for the non-local NN potential U (r, r') 
in quenched QCD at m n ~ 530 MeV. The convergence is tested by the energy and the angular 
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Table 1: The character table of the cubic group O. The notations of classes {E,6C^ 1 ■ ■ ■) and irreducible 
representation (A\,A2, ■ ■ ■) follow Mulliken Symbols. 
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Figure 1: Top: A comparison of the truncated local potential V ( S un (r;£ , ,L = 0) in the 1 Sq channel with PBC 
(blue) and that with APBC (red) at t —t = 9 and r < 2.0 fm. Bottom: A comparison of V^ m (r;E,L) with 
the APBC for the 'So channel (red) and that for the '/^ channel (cyan). 



momentum dependence of the truncated local potential ^""(r^ ,L): To examine the energy de- 
pendence, we have generated VQ Un (r;E,L = 0) in the channel at two different energies, E ~ 
MeV and E ~ 45 MeV. We have found that the energy dependence of the truncated potential is 
very weak. We have also generated VQ Un (r;E,L = 2) for x Di at E ~ 45 MeV in the same setup. 
From the comparison with V^ n {r;E ,L = 0) at E ~ 45 MeV, we have observed that the angular 
momentum dependence of the effective potential is also small. These results indicate that the LO 
potentials constructed from NBS wave functions obtained at E ~ MeV can be safely used for 
< E < 45 MeV and the angular momentum L < 2. 
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An extension of the current analysis to NN potentials in the spin-triplet (5=1) channel is on- 
going. It is an important future task to perform the same analysis for other baryon-baryon potentials 
and for the potentials in full QCD with lighter pion masses. Eventually, we need to extract NLO 
terms such as the spin-orbit potential from the NBS wave functions with finite angular momentum 
on the lattice. 
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